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Abstract 

We study the differential Brancliing ratio and CP asymmetry for the exclusive decay 
B K*l~^l~ in the three Higgs doublet model with additional global 0(2) symmetry in 
the Higgs sector. We analyse dilepton mass square dependency of the these quantities. 
Further, we study the effect of new parameter of the global symmetry in the Higgs sector 
on the differential branching ratio and CP asymmetry. We see that there exist an en- 
hancement in the branching ratio and a considerable CP violation for the relevant process. 
In addition to this, we realize that fixing dilepton mass gives information about the sign of 
the Wilson coefficient C^'^'^. Therefore, the future measurements of the CP asymmetry for 
B — > K*l~^l~ decay will give a powerful information about the sign of Wilson coefficient 
C^^^ and new physics beyond the SM. 
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1 Introduction 



Measurement of the physical quantities for rare B-decays provides sensitive tests for the Stan- 
dard model (SM) and it plays an important role in the determination of the parameters, such 
as Cabbibo-Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants, etc., since 
these decays are induced by flavor changing neutral currents (FCNC) at loop level in the SM. 
Further, they give a comprehensive information in the search of the physics beyond the SM, 
such as, two Higgs Doublet model (2HDM), Minimal Supersymmetric extension of the SM 
(MSSM) 1^, etc. From the experimental point of view, the physical quantities, like Branch- 
ing ratio (Br), CP asymmetry (Acp), forward backward asymmetry (Afb), in rare B-decays 
have an outstanding role to obtain restrictions for the free parameters of the theory under 
consideration. 

Among the rare B-decays, B K*l^l~ , induced by the inclusive process b — > sl^l~, has a 
large Br in the framework of the SM and it can be measured in future experiments. Therefore, 
the study on this process becomes attractive. In the literature, there are various studies on 



these decays in the framework of the SM, 2HDM and MSSM [|[- |T6|. 

CP violating effect is another physical quantity which can give information about the 
free parameters of the model. For B K*l^l^ decay, CP violation almost vanishes in the 
framework of the SM, since the matrix element of the inclusive decay h sl^l~, inducing 
B K*l^l^ process, contains only Vj^VJ*, due to the unitarity of CKM, VibV*, = , {i = u,c,t) 
and smallness of the term VubV*^. This problem was studied in the general 2HDM, so called 



model III, which has a new source for CP violation [|T7|. In that work, the Yukawa couplings are 
taken complex and extra phase angles appear. These new parameters produce a considerable 
CP violation effect for the decay under consideration. 

In this work, we study Br and Acp for the exclusive decay B K*l^l~ in the framework 
of the three Higgs doublet model. Similar to the model III, complex Yukawa couplings are 
possible CP violating sources. However, in the 3HDM, the number of free parameters are 
large compared to that of 2HDM since the Higss sector is extended. We solve this problem by 
introducing a new global 0(2) symmetry in the Higgs sector. 

Even if the theoretical analysis of exclusive decays is more complicated due to the hadronic 
form factors, the experimental investigation of them is easier compared those of inclusive ones. 
Therefore, this work is devoted to the study of the exclusive B K*l^l^ decay. 

The paper is organized as follows: In Section 2, we present our theoretical work based on 
3HDM and the matrix element for the inclusive b — > sl~^l~ {I = e, fi) decay in this model. In 
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Section 3, we calculate Br and Acp for the exclusive decay B — > K*l^l" . Section 4 is devoted 
to discussion and our conclusions. In Appendix, we give some theoretical results for the 3HDM 
and explicit forms of the necessary functions appear in the text. 



The inclusive h 
3HDM 



sl^l decay in the framework of 



In this section, we will derive the matrix element of the inclusive decay h sl~^l~ [l = e,fi), 
which induces the exclusive B K*l^l'' process, in the framework of the 3HDM. We start 
with the general Yukawa interaction. 



Y 



+ pfjQiL<p3UjR + pf^^QiL(l)3DjR + h.c. 



(1) 



where L and R denote chiral projections L{R) = 1/2(1 =i= 75), 0j for i = 1,2, 3, are three scalar 
doublets and rjfj^ , ^fj^, pfj^ are the Yukawa matrices having complex entries, in general. Now, 
we choose scalar Higgs doublets such that the first one describes only the SM part and last two 
carry the information about new physics beyond the SM: 

1 



V2 



, _ 1 / V2H+ 
with the vacuum expectation values. 







V2F^ 



(2) 



< 01 > = 



; < 02 >= ; < 03 >= . 



(3) 



Note that, the similar choice was done in the literature for the general 2HDM ||T8[. The Yukawa 
interaction due to the new physics beyond the SM part is responsible for the Flavor Changing 
(FC) interactions and it can be written as 

jOy,FC = &iL^2UjR + iF,Q^L(|>2DJR + P^QiL^^UjR + pf^QiLfj^sD^R + h.C. . (4) 

Here, the couplings ^^'^ and p^'^ for the charged FC interactions are 



= ^nVc 



p1 



CKM 5 
VcKM C,N , 

Pn Vckm , 

VcKM pN , 



(5) 



and 



(■U,D _ f-,/U,D-.-lcU,DirU,D 
SAT — It^L J S ^ 

pr = (vf'^)-v^'^v^r , (6) 

where the index " N" in denotes the word "neutral". 

At this stage, we obtain the effective Hamiltonian for the inclusive process b —>■ sl^l~ 
by matching the full theory with the effective low energy one at the high scale /i. In this 
calculation, there are additional charged Higgs effects coming from the new charged Higgs 
particles (see eqs. and @). Fortunately, — quark — quark interaction is the same as 
— quark — quark one except new Yukawa couplings and they give additional contributions 
to the Wilson coefficients withouth changing the operator basis Oi{^) (see [0). The Wilson 
coefficients are evaluated from fi down to the lower scale fi ~ 0{m}j) using the renormalization 
group equations. Here the problem is to choose the high scale. In the literature, this scale 
is taken as the mass of charged Higgs, /i = mij±, in the 2HDM, since the evaluation from 
/i = mH± to /i = m\Y gives negligible contribution to the Wilson coefficients( see |]l9l)- In our 
case, there is a new charged Higgs and its mass mp± can be greater compared to mH±. 
However, by introducing a new symmetry in the Higgs sector, we can take that masses of 
and if ^ are the same. Before starting with this discussion, we would like to present the effective 
Hamiltonian, obtained by integrating out the heavy degrees of freedom, here, t quark, W^, H^, 
F^, H^, H"^, and H"^ where H^, F^ and H^,H'^, H^^H^ denote charged and neutral Higgs 
bosons respectively: 

Heff = -i^Vt,VZ i:(a(/i)0.(/i) + CliMifi)) . (7) 

In this equation, are current-current {i = 1, 2, 11, 12), penguin {i = 1, ...6), magnetic penguin 
{i = 7, 8) and semileptonic [i = 9, 10) operators [|T6l , ^ and primed counterparts are their 
flipped chirality partners ||T6|. Ci(/i) and C-{fi) are Wilson coefficients renormalized at the scale 
fi. The initial values of the Wilson coefficients in the SM model, C^'''^'^ {mw), can be found in 
Appendix A. The additional contributions to the initial values of the Wilson coefficients, due 
to two new Higgs scalars are denoted by C^{mw) and for unprimed set of operators we have 

Cl^..6,ll,12("^W') = > 

111'+ Vts ^tb 



mtnib Vis Vtb 



^2 ipN,tt + P*N,tc T/* ) (Pw",tt + P^,tc T/ , 
lih ^ts ^tb 



mtiTib Vts Vtb 



Vtb 

1 T/* T/ 



^ ^ ^^N\tt + CN,tc TT^) (^N,bb + ^N,sbTT-)G2{y) 
mtiii'b <'ts Vtb 

'"t ►^ts ^^t6 

+ -V ^Ptu + PN.tM) iPN,bb + plsb^)G2{y') , 

ribtiibb Vtb 



?7it Vt^ ' Vtb 

+ A + PN,tM) iPN,tt + P%c^)H,{y') , 



TTT,t Vts vtb 



^t Vts Vtb 

and for primed set of operators, 



G'i^...%,ll,12{mw) — , 

''H *^ts vtb 

mtmb Vts Vtb 

'"t Vts Vtb 

I 1 / -*D ^tb N / -(7 I -C/ ^cb \Tp/r\ 

mtUlb Vts Vtb 

C'"{mw) = ^ (G^bsTT? + Cr^ss) i^N,bb + CN,sb^)Giiy) 
'"t Vts Vtb 

+ ZTZT (^N^bs 777 + ^n!ss ) (C^,tt + CN,tc TT ) ^2 (y) 
mtrrib Vts Vtb 



~~2 iP*N\bsT7^ Pn\ss) iPN,bb + PN,sbTr , 

mt Vts Vtb 



+ ^ ^Pt^^V^ + Ptss) iPltt + p1*c17^)G2(i/') , 
mtrrib v4s nfe 



~~2 iPN\bsT7^ + P?^ss) (P^,66 + PN,sbTr , 
^t Vts Vtb 
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C'loimw) - — {Cj^bsT^ + ^N,ss) i^N,bb + ^N,sb7T-)Ll{y) 

^ts ^tb 



+ — iPNMJ^ + Pn,ss) iPN,bb + PN,st^)Lliy') , (9) 

mf ' Vts ' Vtb 
where y = ml / 'mjj± and y' = / m|,± . In eqs. (|) and (j^) we used the redefinition 



= e^'^ . (10) 

The exphcit forms of the functions Fi(2){y), G'i(2)(y), Hiiy) and Li{y) can be found in Appedix 
A. In the calculations, we neglect the contributions of the neutral Higgs bosons to the Wilson 
coefficient C^^^ (see [^). Note that, the neutral Higgs bosons coming from 03 give contribution 



to Cj^^ , including the Yukawa couplings p^^^j and p^^^ {j = d,s ; i = d,s,b), similar to the 



ones coming from 02 [22 



In the 3HDM model, the Higgs sector is extended and this leads to an increase in the number 
of free parameters, namely, masses of new charged and neutral Higgs particles, new Yukawa 
couplings. In the Appendix B, we give the general gauge and CP invariant Higgs potential for 
the 3HDM and present the masses of charged and neutral Higgs particles. Now, our aim is to 
decrease the number of free parameters in the model under consideration. We consider three 
Higgs scalars as orthogonal vectors in a new space, which we call Higgs fiavor space and we 
denote the Higgs fiavor index by "m" , where m = 1, 2, 3. We introduce a new global symmetry 
on the Higgs sector which keeps the 3HDM Lagrangian invariant. Let us take the following 
0(2) transformation: 



2 = COS a 02 + sin a 03 , 

3 = —sin a 02 + cos a 03 , (11) 



where a is the global parameter, which represents a rotation of the vectors 02 and 03 along the 
axis that 0i lies, in the Higgs fiavor space. The kinetic term of the Lagrangian (see Appendix 
B) is invariant under this transformation. The invariance of the potential term can be obtained 
if the following conditions on the free parameters (eq. ^8]) are satisfied: 

C5 = C6 , Cs = Cg , Cn = C12 , 

C2 = C3 = C7 = Clo = , (12) 

and we get 

^(01,02,03) = Ci(0+0i-t;V2)2 + C4[(0+0l- ^72)' + 0^02 + 0^03]' 
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+ C5[(0+0l)(0+02 + 0^03) - (0^02)(0^0l) - (0|03)(0^0l)] 

+ C8([i?e(0+02)]' + [i?e(0+03)]') + Cn([/m(0+02)]' + [/m(0+03)]') 

+ Ci3[/m(0+03)]' + Cl4 (13) 



Therefore, the masses of new particles are 

mp± = mH± = C5 — 
rriH^ = rriH^ = Cg 

V 



2 

2 

V 

niHi = = Cii" 



(14) 



It is the first gain in decreasing the number of free parameters. Now, we apply this transfor- 
mation to the Yukawa Lagrangian (eq.(|l])). This term is invariant if the transformed Yukawa 
matrices satisfy the expressions 

pUiD) fU(D) , -U{D) ■ 

-'U(D) FUiD) ■ . -U{D) /,p.N 

Pij = 'sina + Pi> 'cos a. (15) 

and therefore 

which permits us to parametrize the Yukawa matrices and p^(^) as 

= sin 9 , 

te^stnO , (17) 



P 



where e'^^^^ are real matrices satisfy the equation 

Here T denotes transpose operation. In eq. ([T7|) , we take p^ complex to carry all CP violating 
effects on the third Higgs scalar. 

Finally, we could reduce the number of free parameters, here the Yukawa matrices ^^'^^^ 
and p'^'^^\ by connecting them by the expression given in eq. Further, we take into 

account only the Yukawa couphngs ^'^u^ ^Nbb^ 

and p^bb^ since we assume that the others 



are small due to the discussion given in |16|. Now, we rewrite the contributions of the charged 



Higgs particles to the initial values of the Wilson coefficients as: 
Cl^. . 6,11,12 ("^H/) = , 

+ ^ ^ (^*N,tt ^N,bb + P*N,tt PN,bb) Fiiy) , 

mtmb 

C^imw) = -\{CN%^N,u + PtuPN,u)Giiy) 

+ ZTZr i^*N,ttiN,bb + P*N,ttPN,bb)G2{y) ■, 

rritmi) ' ' ' ' 

C^imw) = -^{CN%^% + pZtPN,u)Hi{y) , 

CUmw) = \{ittlN,u + PtuP\u)Li{y) , (19) 

and neglect the primed coefficients since they include small Yukawa couplings. Here, ^^^^J^b) 
and p^j}ft(pb) ^^"^ obtained by using eq. (0). Note that, with the replacements ^N\tt^N,tt + 
pN,ttPN,tt tN,tt^N,tt and CTy,tt^N,bb + P*N,ttPN,bb tN,tt^N,bb wc get the resuhs for the general 
2HDM (model III) ||16|. By neglecting the primed ones , the initial values of the Wilson 
coefficients can be written as 

Cf^'^^'imw) = Cr\mw)+C!'imw) , (20) 

and using these initial values, we can calculate the coefficients Cf^^^{^) at any lower scale with 
five quark effective theory, namely u, c, d, s, b. In the process under consideration the Wilson 
coefficients Cj^^{fi), Cg^^{fi) and Cio(yu) play the important role in the physical quantities and 
the others enter into expressions with operator mixing. Besides the perturbative part, there 
exist the long distance (LD) effects due to the real cc in the intermediate states, i.e. the cascade 
process B K*il)i K*l^l~ where ^ = 1, ..,6. These effects appear in the Wilson coefficient 
Cg^^ and using a Breit-Wigner form of the resonance propagator |13|, they are added to 
the perturbative one coming from the cc loop: 

C'9"{f^)=Cr\p)+yreson, (21) 

where Yreson in NDR scheme is defined as 

^ 3_ ^ T^T{Vi ll)mv^ 

(3Ci(/i) + C2(/i) + + C4(/i) + 3C5(/i) + C^{^J)) . (22) 
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In eqs. (22), the phenomenological parameter k is taken as k = 2.3 [IC]. The exphcit forms 
of the perturbative parts of the Wilson coefficients Ci{fi) ,i = 1,...,9 including NLO QCD 
corrections can be found in the literature |1T6| , p4t p5[] . 

Finally, neglecting the strange quark mass and primed coefficients, the matrix element for 
b si~^i~ decay is obtained as: 



M 



VM C'/^ifi) S7^(l - 75)& m + CM s^^il - 75)& ilV 



2v^7r 



3 The exclusive B K*l^l decay 

In this section, we study the Branching ratio (Br) and the CP asymmetry {Acp) of the exclusive 
decay B K*l^l~ in the 3HDM. Using the results for the matrix elements {K* |s7^(l ± 75)^! -B), 
and {K* |szo"^i,g^(l ± 75)6! 5) the hadronic matrix element of the B —>■ K*l^l^ decay is 
obtained as [ P7[ : 



M 



2Ae^,p,e>^*g" + iB^e* - iB2{e*q){pB + PK*)t. - iB^{e*q)q^ 



(23) 



where e*^ is the polarization vector of K* meson, pb and Pk* are four momentum vectors of 
B and K* mesons, q = Pb — Pk* and A, C, Bi, and Di i = 1,2,3 are the form factors of the 
relevant process. Their explicit forms can be found in the Appendix C. 
Using eq.(^), we get the double differential decay rate: 

dT G'al^\Vt,V;f X'/' 



dq^dz 



2^'^TT^mB 



2Xm% 



mj,s{l + z')(\A\' + \C\ 



B 



2r 
1 

2r 



m 



Xmlil-z') (iB^f + W) 
'3 {Xil - z'') + 8rs} {\B^f + \D 



- 2Xm%{l - r - s)(l - z"^) {Re {B^BD + Re (DiD;)} 



WbsX^/^z 



{Re {BiC*) + Re {ADD} 



(24) 



where z = cosf3 , (3 is the angle between the momentum of C. lepton and that of B meson in the 
center of mass frame of the lepton pair, A = 1 + + — 2r — 2s — 2rs, r = — ^ and s = — 



rriz 



Acp is another important physical quantity which almost does not exist for the given process 
in the framework of the SM. However, with the choice of the complex Yukawa couplings, it is 



possible that such asymmetry exists, in extended models like 2HDM [|T7|. In our case, the 
model under consideration is the 3HDM with global 0(2) symmetry in the Higgs sector and 
the possible source of CP violation comes from complex Yukawa couplings in the third Higgs 
doublet. Using the definition of Acp 

dr{Bs^K*e+e-) _ dV{Bs^K*e+e-) 
dq^ dq^ 



'^^ dT{Bs-'K*e+e-) . dT(Bs^K'e+e-) " ^ > 

dq^ dq^ 



we get 



i?e(A2)[-2(Pi(/i) + 2P2(/i)) Re{Cl^'{fi)) A + Q 



- -2MA2) _. • (^6) 



In eq. ( PBD we use the same parametrization as in [|T^ 



C',ff{^,) = PM\, + PM , (27) 

where A2 is 

A2 = tt^N bbicos" 9 + %s%r? 6) (28) 

mtrrib ' 

Here the functions Pi(/i) and P2(/^) can be written as the combinations of LO and NLO part, 
namely, 

P,i^^) = P^^(/i) + A^^°(/i), 

P2i^^) = P2^^(/i) + P2^^°(/x), (29) 

and 



P^^(/.) = r^^^/^^P2(y) + ^(V^/^^-V^/^^)G2(i/) 



P^^(/.) = v''^''[C^''irn^) + ^FM] 



zU 12 
2 

mf 

zU 12 



+ QdiC^'^if^) + iVcQ^^(/i)) + g«(— OiT(/i) + iVe— Cf,^(//)) 

777,5 '^fe 

8 

+ C2{mw)T.hiV"' , (30) 

i=l 
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where rj = -^^f^^, hi and a, are numbers appear during the evaluation ||T3|. Pf^'^{fi) is the 
coefficient of A2 in the expression ^^^^Cy^^^^^^ (fi) and P.^^^{fi) is obtained by setting A2 = 
in the same expression. The functions A and Q are defined as 



A = , ^. \A2\i-l -3r + s)+ + v^)'(A - 12r(r - 1)) 



3g2r(l + ^/r) 



(31) 



+ 



8(P,(/.) + P2(/i))P2(/.){^^^T^ + ^r|[A(-4r + .) + 12r(r - 1)^] . 



is 



3mir(-l+r)2 ^ 



+ 2 ( \ ^ s{l-s + 3r)T2Ts\ 32 

Smsq r[—l +r) J 

Here the form factors Ai, A2, Ti, T2, T3 and ^ can be found in the Appendix C. 



4 Discussion 

In the general 3HDM model, there are many free parameters, such as masses of charged and 
neutral Higgs bosons, complex Yukawa couplings, ,pfj^ where i,j are quark fiavor indices. 
The additional global 0(2) symmetry in the Higgs fiavor space forces that the masses of new 
charged Higgs particles to be the same. Further, the masses of the new neutral Higgs particles 
in the second doublet 02, are the same as those of the corresponding ones in the third doublet 
03. This symmetry also connects the Yukawa matrices in the second and third doublet (see eq. 
(p!7|)). The Yukawa couplings, which are entries of Yukawa matrices, can be restricted using 
the experimental measurements. In our calculations, we neglect all Yukawa couplings except 
^^j^ and ^§1,1, tt ^"^^ pNbb respecting the CLEO measurement announced recently PS| , 

Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10~^ , (33) 

This section is devoted to the study of the dependencies of Br and Acp for the decay 
B — > K*l^l~ , for the selected parameters of the 3HDM {e^ui ^nw ^^e angle 9) with 
0(2) symmetry in the Higgs sector. In our analysis, we restricted lOy'^'^^l in the region 0.257 < 
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\C^^^ < 0.439, coming from CLEO measurement where upper and lower limits were 



calculated in following the procedure given in p9[. This restriction allows us to define a 
constraint region for the parameter e^^^ in terms of e^^f, and 6. Our numerical calculations 
based on this restriction and the constraint for the angle 6 due to the experimental upper 
limit of neutron electric dipole moment, namely dn < 10~^^e-cm which places a upper bound 



on the couplings with the expression: 



^Ntt^Nhb)sin^d < 1-0 for Mh± ~ 200 GeV ^ 



mtmt ^^N,tt <^Af,66 

Throughout these calculations, we take the charged Higgs mass mF± 
scale /i = rrifo and we use the input values given in Table (0). 



mH± = 400 GeV, the 



Parameter 


Value 


rric 


1.4 (GeV) 


nib 


4.8 (GeV) 


em 


129 




0.04 


rriB^ 


5.28 (GeV) 
3.9610-3 


^tot{Bd) 




175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


^QCD 


0.214 (GeV) 


as{niz) 


0.117 


sinOw 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 



In fig. |l], we plot the differential Br of the decay B K*l^l with respect to the dilepton 
mass for e^bb — 40 m^, sinO = 0.5 and charged Higgs mass nijj± = 400 GeV at the scale 
fi = nib- This figure represents the case where the ratio \rtb\ = I^tt^I < 1. Here the differential 



Br lies in the region bounded by solid lines for C^-^-^ > and by dashed lines for C^^^ < 0. It 
is shown that there is an enhancement for Cj-^'^ > case compared to the SM (dotted line). 
Further, the restriction region of the differential Br for C^^^ > case is broader than the one 
for Cj^^ < 0. Fig. 1^ is devoted the dependence of the differential Br to sin 9 for \rtb\ < 1, 
^Nbb = 40mfe, = 12Gel^^ and charged Higgs mass niH± = 400 Gel^ at the scale fi = nib- 
Here, the differential Br lies in the region between solid lines for Cj^^ > and lies in the region 
between dashed lines for Cj^^ < 0. For Cj-^-^ > 0, there is a weak dependence to sin 9 especially 
for sin 9 < 0.5. For Cy-^'^ < this dependence almost vanishes. Now, we present the Br for 
three different phase angles {sin9 = 0.1, 0, 5, 0.9) in two different dilepton mass regions (Table 



N,bb 
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sinO 


Cf ^ > 


Cj^^ < 


2 

Q regions 


fl 1 


1 80 1 0"'' < Br < 2 7^ ^ 0^^ 


95 1 O"*' < Rr < 1 07 1 0^^ 


I 




0.9610"'^ < 5r < 1.07 lO"*^ 


0.6610"** <Br< 0.7210"** 


II 


0.5 


1.7210-'^ < 5r < 2.1210-*^ 


0.9410"** <Br< 1.0610"** 


I 




0.9410-^ <Br< 1.0510-*^ 


0.6610"** <Br< 0.7210"** 


II 


0.9 


1.0410-** <Br< 1.26 10"^ 


1.0410"** <Br< 1.0510"** 


I 




0.71 10-** <Br< 0.77 10"^ 


0.6910"** <Br< 0.7110"** 


II 



Table 2: Br for regions I ( 1 GeV < < mj/^ - 20 MeV ) and II (mj/^ + 20 MeV < < 
m^/ - 20 MeV ) 

In figs. D (^) we plot Acp of the decay B K*l^l~ with respect to the dilepton mass 
square, g^, for e^^^f, = 40 m^, sinO = 0.1 {sin9 = 0.5) in the case where the ratio \rth\ < 1- For 
Cj^^ > 0, Acp is restricted in the region bounded by solid lines and for Cy-^-^ < it lies between 
dashed lines. Acp changes sign almost at the value g^ ~ 9 GeV^ for C^^^ > case. However, 
for Cj^^ < 0, it can have both signs for any g^ value. Acp enhances strongly with increasing 
value of sin 6. For completeness, we present the average value of Acp for three different phase 
angles {sin6 = 0.1, 0, 5, 0.9) in two different dilepton mass regions (Table |^), 



sinO 


C'^^^ > 


C," < 


q'^ regions 


0.1 


-3.35 10""* < Acp < -0.74 10"* 


-0.74 10"* < Acp < 1.71 10""^ 


I 




0.43 10""* < Acp < 3.2210"* 


-1.05 10"* < Acp < 0.43 10"* 


II 


0.5 


-1.0510"^ < Acp < -1.0410"^ 


-0.25 10"^ < Acp < 0.55 10"^ 


I 




0.9710"^ < Acp < 1.1710"^ 


-0.3410"^ < Acp < 0.1410"^ 


II 


0.9 


-3.32 10"^ < Acp < -0.91 10"^ 


-0.9010"^ < Acp < 2.7210"^ 


I 




0.5310"^ < Acp < 2.6610"^ 


-1.8610"^ < Acp < 0.5310"^ 


II 



Table 3: The average CP asymmetry Acp for regions I ( 1 GeV < y/f < mj/^ - 20 MeV ) 
and II (m j/^ + 20 MeV <^< m^, - 20 MeV ) 

Figs. 1^ and ^ are devoted to sinQ dependence of Acp for g^ = 'IGeV'^ and g^ = 12GeV'^ 
respectively. Here, Acp lies in the region bounded by solid lines for G^-^-^ > or by dashed 
lines for Gj^-^ < 0. With decreasing sin 6, Acp decreases as expected and the restriction region 
becomes narrower, for both Gj'^^ > and Gj^^ < 0. Further, for fixed g^ values, the sign of 
Acp does not change with changing sin 9 for Gj'^^ > , contrary to the case G^-^-^ < 0. This 
is informative in the determination of the sign of Gj^-^ with the experimental measurement of 
Acp at fixed g^. Note that the similar situation exist for the general 2HDM with complex 
Yukawa couplings (see |T7[|). 
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Now we would like to summarize our results: 

• Br for the process under consideration is at the order of 10^® for \rtb\ < 1 and it is greater 
for C^-^-^ > compared to C^-^^ < 0. Further, it is not sensitive to sin 6 especially for 
C'^" < 0. 

• \Acp\ increases with increasing sin6. For C^-^-^ > 0, Aqp changes sign at the value, 

~ 9GeV^, however it can have any sign for Cj-^-^ < 0. For the case \rtb\ » 1, Acp 
almost vanishes (~ 10^^^) since sin 6 should be small due to the restriction coming from 
the limit on neutron electric dipole moment. 

• For the fixed value of q"^ and Cj'^-^ > 0, Acp can be either negative or positive when sin9 
varies. For C^^^ < 0, it can have both signs. This shows that with the measurement of 
Acp for fixed q^, it is possible to detect the sign of C^^^ , which is an interesting result. 

Therefore, the experimental investigation of Acp ensure a crucial test for new physics and 
also the sign of C^^^^ ■ 
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Appendix 



A The Wilson coefficients in the SM and the functions 
appear in these coefficients 

The initial values of the Wilson coefficients for the relevant process in the SM are 



Cl,3,...6,ll,12("^W/) 

and the primed ones are 



, 

1 , 



4(s - 1)4 

4(x - ly 
1 



In X H ; r:^ 

2A{x - 1)'^ 
—x^ + bx^ + 2x 



In X + 



siiT?9w 
1 



B{x) + 



8(x - 1)3 
4 sin^ Ow 



w 



C{x)-D{x) + -, , 



sin^ 6\Y 



{B{x)-C{x)) 



The functions appear in these coefficients are 



B{x) 
C{x) 
D{x) 



-X 



+ 



X 



Inx 



X — 1 (x — 1)^ 
x/2 - 3 ^ 3x/2 + 1 



Inx 



X — 1 ' (x — 1)^ 
19x736 + 25x736 -x76 + 5x73 - 3x2 ^ ^g^/g _ 4/g 



(x - 1)3 

and in the coefficients C^'^^ (eqs. (|) and 
y{7-5y-8y') , y\3y - 2) 



{x-iy 



Inx 



are 



Fi{y) 
F2{y) 



72(?/-l)3 +12(?/-1)4 
2/(51/ -3) , y{-3y + 2) 



hiy , 



12(2/ -1)2 6(y-l)3 
2/(-?/2 + 5?/ + 2) , -1 



Iny , 



24(y - 1)= 



4(2/ - 1)^ 



\ny , 



(34) 



(35) 



(36) 
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G2{y) 
Hi{y) 
Li{y) 



yjy - 3) y_ 

4{y - 1)2 1{y - 
1 — ALsirx^dw xy 

sin^Ow 8 

1 xy 



1)3 
1 



y 



siri^Ow 8 



+ 



(1/ 
1 



y - 1 (y - 1)' 



1)2 

Iny 



hiy 



47|/2 - 79|/ + 38 3|/3 - 
108(y- 1)3 1% 



6y + 4 



(37) 



B 3 Higgs doublet model 

We consider three complex, SU{2) doublet scalar fields 0i (i = 1,2,3). The gauge and CP 
invariant Higgs potential which spontaneously breaks SU{2) x U{1) down to U{1) can be 
written as: 

y(01,02,03) = Ci(0+0i-^;V2)2 + C2(0+02)' 

+ C3(03+03)' + C4[(0|01 - ^V2) + 0^02 + <Pt<p3? 
+ C5[(0M)(0^02)-(0^02)(02+0l)] 
+ C6[(0M)(03+03)-(0l^03)(03'"0l)] 
+ C7[(0+02)(0^03)-(02V3)(03'"02)] 

+ C8[i?e(0+02)]' + C9[i?e(0+03)]' + Cio[i?e(02+03)]' 

+ cii[/m(0+02)]' + ci2[/m(0+03)]' + ci3[/m(0+03)]2 + ci4 (38) 

Here, we assume that only 0i has vacuum expectation value (see section 2). The parameters 
Ci arc real to ensure the hcrmiticity of the potential term. Further, the Higgs sector does 
not violate CP and all possible CP violation effects are based on the choice of the Yukawa 
couplings. 

The Higgs boson squared mass matrix can be obtained by the expression 



niij = 



d0. de, 



(39) 



Here 6i are real fields satisfying 



^1 + 192 , 

O7, 
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= ^12, (40) 

where x'^^ X° ^ind are the SM particles and H^, F^, W {i = 1,2,3,4) are new particles 
existing in 3HDM (see eq. (|^) ). Note that, these fields are mass eigenstates, thanks to choice 
eq.(^). Diagonalizing this matrix, we get masses of new Higgs particles as: 







mH± = 


2 ^ 










mF± = 


2 ^ 










mm = 












mm = 












mm = 


— Co 

2 ^ 










mm = 









(41) 

In eq. (^OD, X° are golstone bosons, which can be eaten up in unitary gauge, is the 

SM Higgs which has mass mm = 2v^{ci + C4). H^, are scalar and H^, are pseudoscalar 
particles due to new physics. Note that and are denoted by and in the literature. 
For completeness, we also present the kinetic term for the 3HDM: 

{dyi-t^-By,-t^-<P,^Wn (42) 

where 

0^= ( ^0 ) ^ = 1,2,3. (43) 

C The form factors for the decay B K*l^l~ 

The structure functions appear in eq. ( p3| ) are 

A = -d^^ 4C,^^^^Ti, 

mB + mx* q 
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-Claims + m^O^i - ^C'r^'^iml - ml.)T, 



B2 = 




B3 = 




c = 


n 




Di = 




D2 = 


n 


D3 = 


C ^ 

'-'10 



A, 



+ rriK 



AG 



eff^b 
7 ~ 



T2 + 



2 2-^3 



V 



A. 



-{A,- A,) . 



(44) 



We use the dependent expression which is calculated in the framework of light-cone QCD 
sum rules in to calculate the hadronic formfactors V, Ai, A2, Aq, Ti, T2 and T3: 

m 



(45) 



where the values of parameters F{0), ap and hp are listed in Table 4. 





m 


ap 


bp 


Ai 


0.34 ±0.05 


0.60 


-0.023 


A2 


0.28 ±0.04 


1.18 


0.281 


V 


0.46 ±0.07 


1.55 


0.575 


Ti 


0.19 ±0.03 


1.59 


0.615 


T2 


0.19 ±0.03 


0.49 


-0.241 


T3 


0.13 ±0.02 


1.20 


0.098 



Table 4: The values of parameters existing in eq. (|45|) for the various form factors of the tran- 
sition B K*. 
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Figure 1: Differential Br as a function of for sin 9 = 0.5, 6^55 = 40 in tlie region 
\rtb\ < 1, at the scale /i = mt, including LD effects. Here differential Br is restricted in the 
region bounded by solid lines for C^^-^ > and by dashed lines for C^-^^ < 0. Dotted line 
represents the SM result withouth LD effects. 
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Figure 2: Differential Br as a function of sin 6 for q"^ = 12GeV'^, (^^bb — 40 in the region 
\rtb\ < 1, at the scale = rrih, including LD effects. Here differential Br is restricted in the 
region bounded by solid lines for Cy'^'^ > and by dashed lines for C^^^ < 0. Dotted line 
represents the SM result withouth LD effects. 
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Figure 3: Acp as a function of for sinO — 0.1, e^;,^ = 40 in the region \rtb\ < 1, at the 
scale fi = mil, including LD effects. Here Acp is restricted in the region bounded by solid lines 
for C^'^'^ > and by dashed lines for Cy'^'^ < . 
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Figure 4: The same as Fig|^, but for sin 9 = 0.5. 
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Figure 5: Aqp as a function of sin 9 for = AGeV"^, e^^^ = 40 m;, in the region \rtb\ < 1, at 
the scale /i = m^, including LD effects. Here Acp is restricted in the region bounded by solid 
lines for C^^^ > and by dashed lines for C^^^ < 0. 
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